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Abstract. We compute the rates P at which acoustic energy is injected into the solar radial p modes for several 
solar models. The solar models are compu ted with two differe nt local treatments of convection: the classical 
mixing-length theory (MLT hereafter) and ICanuto et all il996l CGM96 hereafter) 's formulation. Among the 
models investigated here, our best models reproduce both the solar radius and the solar luminosity at solar 
age and the observed Balmer line profiles. For the MLT treatment, the rates P do depend significantly on the 
properties of the atmosphere whereas for the CGM96's treatment the dependence of P on the properties of the 
atmosphere is found smaller than the error bars attached to the seismic measurements. The excitation rates P 
for modes associated with the MLT models are significantly underestimated compared with the solar seismic 
constraints. The CGM models yield values for P closer to the seismic data than the MLT models. We conclude 
that the solar p-mode excitation rates provide valuable constraints and according to the present investigation 
clearly favor the CGM treatment with respect to the MLT, although neither of them yields values of P as close 
to the observations as recently found for 3D numerical simulations. 
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1. Introduction 

In the outermost part of the convective zone (CZ) of in- 
termediate mass stars, convection is highly superadiabatic 
because of the rapid radiative heat gains and losses of the 
convective fluid. In that region, entropy fluctuations are 
the largest and the resulting decrease of the convective 
transport efficiency is compensated by a large increase of 
the eddy motions which is responsible for the oscillation 
mode driving. Modelling inefficient convection is complex. 
3D numerical simulations are now being performed but re- 
main still very time consuming. Hence for massive stellar 
computations, ID stellar models are used in which only 
simple prescriptions of convection are implemented. 

Among the se simplified treatments, 

ICanuto fc Mazzitellil l|l99l[ CM91 hereafter) 's ap- 
proach differs from the classical mixing length approach 
(MLT hereafter) in that it takes into account the con- 
tribution of eddies with different sizes in the calculation 
of the convective flux and velocity while keeping the 
computational expenses a s low as the M L T. An improved 
version was proposed by ICanuto et alJ l)l996l CGM96 



hereafter) which takes into account the feedback of 
the turbulence on the energy input from the source 
which generates turbulent convection. These multi-eddy 
convection models are usually refered to as Full Spectrum 
of Turbulence (FST) models. 

Several non-local formulat ions of convection have also 
been propo sed lGough.l977t IXiondll978l Il985t ICanutd 
1992, I99J; ICanuto fc Dubovikovlll998(l . However, we fo- 
cus here on the effects of proposed improvements in the 
description of the energy spectrum and therefore consider 
only local treatments and compare FST models with MLT 
ones. 

Any model of convection must satisfy several 
observational constraints provided by our Sun: the 
solar radius at the solar age, the Balmer line profiles 
and the uvby color indices. The MLT, CM91's and 
CGM96's local treatments have bee n confronted to 
these qbsery at ional constraints (e.g. iFuhrmann et alJ 
199^ IT99i fvan't Veer-Menneret fc Megessieij Il99fit 
Smallev fc Kupkal h()^ iFiornkopf. [iM Iw^it^r et a,lJ 
2OO2I: iMontalban et alJ l2004l) . One main result is that 
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these observational quantities are more sensitive to the 
adopted value of the convective scale length of the eddies 
than to the formulation of convection. 
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Solar seismic observations provide strong additional 
constraints. Comparisons of theoretical oscillation fre- 
quencies with observed solar ones have shown for instance 
that significant improvement in the agreement between 
observation and model at high frequen cy and degree £ 
can b e achieved with 3D simulations I Rosenthal et al.l 
We are here interested in amplitudes of oscillation 
which can also bring several constraints on the convective 
process in the outer solar envelope. Indeed, the ampli- 
tudes of solar-like oscillations result from a balance be- 
tween excitation and damping. Measurements of the os- 
cillation mode growth rates (through their line-widths) 
and of the mode amplitudes enable the evaluation of 
the excitation rates P. Excitation of solar-like oscilla- 
tions is known to be both due to turbulent convective 
motions thro ugh the driving by the turbulent Reynolds 
stresses (sec 'Coldreich fc KeelevI Il977t iBalmforthI Il992t 
ISamadi fc Goupi l .2001) and due to the advection of tur- 
bule nt entropy fluctuations by the turbulent movements 
(see ISamadi fc Goupill l200lh . The excitation rates P are 
thus directly related to the velocity of the convective ele- 
ments and to the amount of thermal energy advected by 
convective motions (i.e. the convective flux). The excita- 
tion rates then crucially depend on t he way the convectiv e 
velocity and flux are modelled (see iHoudek et al]ll999^ . 
Solar seismic measurements therefore provide - through a 
model of mode excitation - additional constraints on the 
stellar convective properties. In this framework, the goal 
of the present paper is to investigate the influence of dif- 
ferent local treatments of convection on the calculation of 
the rates at which energy is injected into the solar radial 
p modes and to compare our results with the solar seismic 
constraints. 

For this purpose we compute two calibrated solar 
models with the Bohm-Vitense formulation of the MLT 
(Bohm-Vitcnse 1958, hereafter BV) and with the CGM96 
multi-eddy convection treatment. In each case, the same 
convection formulation is adopted for the interior and 
the model atmosphere. Models for the internal struc- 
ture are built so as to reproduce the solar radius and 
the solar luminosity at the solar age. The atmosphere 
of each model is constructed using a T(t) la w whic h 
is de rived from a Kurucz's model atmosphere |Kiiruczl 
I1993I) computed with the same convection formulation 
(as described in Heiter ct al. 2002) . These model atmo- 
spheres are built in order to provide the best agreement 
between synthetic and observed Balmer line profiles (as 
in fvan't V cer-Menncrct & Mcgessier 1996, for the MLT 
treatment) (Sect. 2). The matching of the model atmo- 
sphere with the interior model is performed - in the man- 
ner of |Morel et al. ( 1994) - by ensuring the continuity of 
the temperature gradient, V, and of the convective flux 
in a transition region between the interior and the atmo- 
sphere. 

We compute also two models with an Eddington gray 
atmosphere. One with the MLT treatment and the sec- 
ond one with the CGM96 formulation. These two addi- 
tional models are considered for comparison purpose only. 



Indeed, they have an atmosphere with the same mixing- 
length parameter as in the interior and do not reproduce 
the Balmer line profiles. As a consequence, in contrast 
with the interior models including a Kurucz's atmosphere 
as described above, their atmospheres do not fulfill con- 
straints on the properties of the convection at the surface. 

Calculation of the excitation rates requires the com- 
putation of the convective flux, Fc, and of the convective 
velocity, v. This is done in Sect. 3 by paying special atten- 
tion to the problem of the transition region. Indeed, the 
continuity of V and of Fc through the transition region 
imposes a spatial variation of the mixing-length param- 
eter in the transition region. This variable mixing-length 
parameter is then used in Sect. 3 to compute v. Note that 
our approach, which is used here to c o mput e v and P, 
is different from that of ISchlattl et alJ l)l99lt) who built 
stellar models assuming a spatially varying mixing-length 
parameter, with a spatial variation imposed a priori from 
a comparison to 2D numerical simulations of convection, 
in order to compute p mode frequencies. 

As a last step (Sect. 4), we compute the adiabatic 
eigenmodes and the excitation rates P f or each model. The 
adopt ed model of excitation is that of ISamadi fc Goupil 
(j200li Paper I hereafter) in which the characteris- 
tic wavenumber fcg, the wavenumber dependency of the 
turbulent spectra as well as the frequency component 
iXk) of the correlation product of the turbulent veloc- 
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with solar seismic constraints then allows to conclude 
about the best local treatment of convection in the solar 
case (Sect. 5). 

2. Solar models 

All solar mod els discussed here are computed with the 
CESAM code jMorellll997(l including the following input 
physics and numerical features: 

1. Equation of state (EOS): CEFF EOS 

|Christc nsen-Dalsgaard fc Dappenlll992l). 

2. Opacities: OP AL lllglesias fc Roger^ 1 ^961) d ata com- 

data for 

given for 



plemented by [Alexander fc Ferguson (|1 
T < 10^ K, both s et of data being 
iGrevesse fc Noelsl lll993l) solar mixture. 
Thermonuclear reaction rates: ICaughlan fc Fowled 
(1988). 

Convection: either MLT or CGM96's formalism. The 
same convection formalism has been used in the inte- 
rior and in the model atmosphere. 
Microscopic diffusion: all models include micro- 
scopic diffusion of helium and heavy elements cal- 
culated according to th e simplified formalism of 
iMichaiid fc ProffittI ^^1993^ where heavy elements are 
treated as trace elements. 

Chemical composi tion and mixing length parameter 
for convection: the IGrevesse fc Noelsl (I1993I) heavy el- 
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ements solar mixture has been adopted. The constraint 
that solar models have the observed solar luminosity 
and radius at solar age yields the initial helium con- 
tent Yq and the mixing length parameter of the in- 
terior model Qfj (solar model calibration). Microscopic 
diffusion modifies the surface composition, therefore 
the initial ratio of heavy elements to hydrogen (Z/X)o 
is adjusted so as to get the ratio {Z/X)q = 0.0245 at 
solar age. 

7. The models were calculated with 285 shells in the at- 
mosphere and about 2000 shells in the interior. 

The CGM96 fo rmulation o f conv ection is imple- 
mented acco r ding t o lHeiter et al. I ll2002l) . In contrast with 
Heiter et al.l l|2002l) . we use for the two formulations a 
characteristic scale length of convection which is the 
mixing-length A = a Hp where Hp is the pressure scale 
heigth and a is the mixing-length parameter which can 
be different in the interior and in the atmosphere. 

2.1. The "Kurucz models" (KMLT and KCGM 
models) 

We consider here two stellar models: one computed with 
the MLT formulation of convection and the second one 
with the CGM96 formulation. They will be hereafter la- 
belled as KMLT model and KCGM model respectively. 

Treatment of the atmosphere: The model atmospheres 
of those models are computed using the A TLAS 9 code 
ijKu rucz 1992t) described in Heiter et al.j l|2002l) . Solar 
model atmospheres are built assuming different values 
for tta, the mixing-length parameter assumed for the 
model atmosphere: = 0.4,0.5,0.6,0.7. The model at- 
mospheres with Q!a — 0.4 and — 0.5 provide the best 
agreement between synthetic and observed Balmer 
line profiles for the two formulations of convection. This 
is shown in Fig. ^ for the CGM model atmosphere . For 
the MLT tre atment see 'Fuhr mann et alJ lll993lll994l) and 
Ivan't Veer-M cnneret & Megc ssieil ()l996|) . Above aa — 
0.6, the synthetic profile rapidly departs from the ob- 
served one as well as the effective temperature Toff from 
the known solar T^g. 

There are no significant differences for the Balmer 
line profile between the model atmospheres with a^. = 0.4 
and tta = 0.5. Among those model atmospheres we adopt 
arbitrary those with — 0.5. Indeed, choosing the model 
atmospheres with eta — 0.4 instead of oia = 0.5 will not 
change the conclusions of this article. 

For each formulation of convection we then obtain a 
T(t) law. 

The atmospheres of the KMLT and KCGM stellar 
models ar e recomputed accor ding to the procedure de- 
scribed in iMorel et alJ (^^) from the T(r)-laws men- 
tioned above: the fit between interior (where the diffusion 
approximation is valid) and atmosphere is performed in a 
region where ti <, t ^ T2 (an acceptable range of values 
for Ti and T2 is discussed below). In the interior region, 
where r <; r2, the temperature gradient Vi is obtained 



from the MLT or CGM96 formalism. In the atmospheric 
region, where r ^ ri , the temperature gradient Va is com- 
puted using the T — t law of the model atmosphere built 
with the same model of convection as in the interior. In the 
transition region, where ti r ^ T2, in order to ensure 
the continuity of the temperature gradient, V is obtained 
by a linear interpolation of Vi and Va as a function of the 
optical depth as follows: 

V = /3(r)Va + (1 - /3(T))Vi (1) 

where /3(r) = {t2 -t)/{t2-ti). 

Once the temperature gradients of the interior and the 
atmosphere are linked together in the transition region ac- 
cording to Eq. we compute afterward in that region 
the convective flux and an equivalent mixing-length pa- 
rameter (i.e. a depth dependent mixing-length parameter) 
as explained in Sect. 3. 

Acceptable ranges for ti and T2: Using a Newton- 
Raphson scheme, Vi is adjusted in order that -fi + -F'/ad ~ 
LQ/Anr^, where Lq is the luminosity of the Sun and fJ^^j 
is the radiative flux of the internal model. Calculation of 
Fc"^ assumes for Fj-ad the diffusion approximation for the 
radiative transfer. This approximatio n is valid at rathe r 
high values of t, typically r > 10 fsee lMorel et alJll994j) . 
Therefore, ti cannot be much smaller than r ~ 10. 
Otherwise, Vi will have an unrealistic contribution to V 
below T ~ 10. On the other hand, for the calculation of 
F^, the radiative flux Fj-ad is based on a Kurucz 's model 
atmosphere which treats the radiative transfer more real- 
istically than the diffusion approximation. 

The Kurucz model atmosphere is based on iKurucd 
((1992, 1993) opacity tables which are given up to T ~ 
2 • 10^ K and P = 10^ dyn cm~^. As a consequence, 
T2 cannot be larger than r ~ 10^'^, i.e. layers for which 
T ~ 30000 K. In order to ensure a satisfactory continuity 
of the temperature gradient T2 must be sufficiently larger 
than Ti. On the other hand, the transition region should 
be as small as possible, i.e. T2 — ti must be as small as pos- 
sible. The main constraint for this region is thus to avoid 
discontinuities between the interior and the atmosphere. 
It is defined through an empirical procedure rather than 
based on a strict physical theory. 

In practice, we find that ri = 4 is the minimal accept- 
able value for ti; below this value the bias introduced by 
the diffusion approximation has a significant effect on Fc. 
In addition, we find that above T2 ~ 50, the convective 
velocity v (see Sect. I3.2|l shows a pronounced "kink" at 
r = T2 for the CGM model (see Fig. 13 such a "kink" 
is also observed for Fc, but it is less pronounced). On the 
other hand, the choice of T2 = 20 avoids the angular point. 
For the MLT model, whatever the value ofT2,v shows such 
a feature. This is a consequence of the much larger values 
of the mixing length parameter required in the interior 
in order to still obtain the correct solar radius, if a lower 
value of a is used also in regions further within the enve- 
lope. The requirement of matching Rq hence provides a 
more stringent upper limit for the choice of T2 (cf. also the 
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Fig. 1. The observed solar profile is compared to theoretical ones computed with CGM models and assuming 
different values for aa- Abscissae are distances in A from the line center, ordinates are the flux in the profile normalized 
to the continuum. The large scatter is due to the presence of many spectral lines which overlap and cause an apparent 
enlargement of the true profile. 



discussion in lMontalban et al .120041 on the computation of 
solar entropy as a function of radius) . In the following, we 
will consider ri = 4 and T2 = 20 as our optimal choice. 

Calibration: The mixing-length parameter ai for the 
internal structure, Yq, and {Z/X)q are adjusted such 
that the stellar model simultaneously reproduces the so- 
lar radius, the solar luminosity, and the observed ratio 
{Z/X)q = 0.0245 at the solar age. The calibration yields 
{Z/X)o = 0.0279 and Yq = 0.275. At solar age, the he- 
lium abundance in the convection zone is y = 0.246 in 
reasonable agreement with the value Y = 0.249 ± 0.003 
obtained from seismology l)Basulll9 97). Table Ogives the 
calibrated values of the mixing-length parameters ai for 
each complete solar model. The radius resulting from the 
adjustement of ai as well as the size of the convective zone 
are given in Table El All the interior models have a depth 
of the convective zone of ~ 0.286 Rq which is in good 
agreement wit h the value of 0.287 ± 0.003 Rp, determined 
seismically by IChristensen-Dalseaard et all |l99l|l . 



2.2. Eddington approximation based models (ElVILT 
and ECGM models) 

For comparison purpose, we consider here two additional 
stellar models with an Eddington classical gray atmo- 
sphere instead of the Kurucz atmosphere models described 
in Sect l2.1l One of these models assumes the MLT formu- 
lation of convection and the other one the CGM96 for- 
mulation. In the following they will be labelled as EMLT 

Table 1. Values of the mixing-length parameters of the 
KCGM and KMLT models: ai (for the interior) and 
(for the model atmosphere) . ai results from the calibration 
of the full model while ag, is fixed (see Sect. 2.1). 

model ai Qa 
KMLT 2^51 050" 
KCGM 0.78 0.50 
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model and ECGM model, respectively. The mixing-length 
parameter a of these models (the same a in the interior as 
in the atmosphere) is adjusted in order to reproduce the 
solar luminosity and radius at the solar age. However, as 
mentioned in the introduction, these models do not repro- 
duce the Balmer line profiles. Table 13 gives the calibrated 
values of the mixing-length parameters. 

2.3. Comments 

With the CGM96 treatment, ai is found less than one and 
closer to = 0.5. In contrast, with the MLT treatment 
the value of ai is much larger than aa — 0.5. 

The CGM models result in a much smaller value for 
the mixing-length parameter (ai for the KCGM model and 
a for the ECGM model) than the MLT models because 
the convection in nearly adiabatic regions is more efficient 
with the CGM96 formulation than with the MLT one. 
Indeed, for the same value of the mixing-length parame- 
ter and in the region of high convective efficiency (below 
the superadiabatic region), the CGM96 treatment predicts 
a convective flux ten times larger than the MLT one for 
a given superadiabatic gradient. In a solar model this be- 
havior results in a gradient closer to the adiabatic one 
below a sm aller superadiabatic zone in comparison to the 
MLT case ijCanuto et al.lll99fi >). At the top of the quasi- 
adiabatic region, energy is predominantly transported by 
convection, such that Fc — F where Fc is the convec- 
tive flux and F is the total energy flux. Above, convective 
transport is no longer efficient {Fc < F). Therefore, the so- 
lar energy flux F at the top of the quasi-adiabatic region 
can be reproduced for the CGM models with a smaller 
value of the mixing-length parameter than that one re- 
quired for the MLT. 

The superadiabatic gradients are displayed in Fig.[21as 
a function of t. The angular point observed in V — Vad at 



Table 2. AR = R — Rq, where Rq is the radius at the 
photosphere (i.e. at T = Toff), and depth of the convec- 
tive zone (CZ) for the KCGM and KMLT models. These 
quantities are given with respect to the solar radius Rq 
(we assume ICuenther et al ] lll992l) 's value of Rq). 



model AR/Rq depth CZ 



the optical depth r ~ 20 for the KMLT model corresponds 
to the matching point {i.e. at Rq — r ~ 100 km, where 
Rq is the solar radius at the photosphere, defined to be 
where r = 2/3 and calibrated at the precision level given 
in Table ^ . It results from the large difference between 
ai and eta- This difference is much smaller for the KCGM 
model and therefore the "kink" at r ~ 20 is much less 
pronounced. 

Our last comment concerns the large difference in the 
value of the mixing-length parameter between the KMLT 
and the EMLT models: Both models differs only by the 
structure of their upper most layers located 50 - 100km 
below the photosphere (see Fig. 3, middle) which represent 
a tiny fraction of the convection zone depth. The KMLT 
has a model atmosphere in which convection is much less 
efficient than that of the EMLT model as a consequence of 
the fact that aa = 0.5 in the Kurucz's model atmosphere 
(see Sec. 2.1). This is the reason why the superadiabatic 
gradient, V — Vad, is reaching at that depth much higher 
values for the KMLT model than for the EMLT model. 
Now if V — Vad from two models are vastly different, so 
is the entropy jump of both. Hence, if a certain entropy 
of the interior convection zone and thus a certain radius 
of the convection zone shall be matched, a much more 
drastic change in a is required with the KMLT model to 
av oid a too large en t ropy j ump (see a detailled discussion 
in lMontalban et all l)2004 '). 
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Table 3. Values of the the mixing-length parameter a 
of the ECGM and EMLT models obtained for cafibrated 
solar models. 

model a 
EMLT L76 
ECGM 0.69 



Fig. 2. The superadiabatic gradient V — Vad is plotted 
versus the optical depth r in the outer region for the 
KMLT model (dashed line), EMLT model (dot-dashed 
line), ECGM model (dot-dot-dot-dashed line) and KCGM 
model (solid line). The dotted vertical lines correspond to 
the layers where r = ri = 4 and r = r2 = 20, respectively, 
and delimit the transition region. 
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3. Convective velocity and entropy fluctuations 

3.1. Convective flux 

Part of the mode excitation rates stems from the advec- 
tion of turbulent entropy fluctuations by turbulent mo- 
tions (the so-called "entropy source term"). This term 
scales - see Paper I - as < > u§ where st rep- 
resents the entropy fluctuations due to turbulent convec- 
tion, as = {dp/ds)p, s, p and p are respectively the den- 
sity, the entropy and the pressure, <> denotes spatial and 
time average, Uq = 1/3 < >, where u is the velocity 
field associated with the turbulence. The factor 1/3 arises 
from the simplifying assumption made in Paper I that the 
acoustic emission is injected into the p-mode isotropically 
in all three directions. 

The "entropy source term" scales as the square of the 
convective flux Fc - Indeed, we show (see Appendix A) that 
it scales as {as/ poT^f ($/3) F^, where Tq and po are the 
mean temperature and density, respectively, and ^ is the 
spatially aver aged anisotropy factor which is defined as 
l|noughlll97i : 

- ^ - 4 (2) 
<ul> 

where Uz is the vertical component of u and v and w are 
defined as =< v?' > and =< >, respectively. 

For the CGM96 formul ation, Fg is compute d according 
to Eqs. (2) and (17)-(21) in lHeiter et"al] l(2002l) and for the 
M LT treatme nt, it is calculated according to Eq. (14.118) 
inlc^llMi Chap. 14). 

Fc can be viewed as function of a and V, Fc = /i(V, a), 
where h is given by the adopted model of convection. 

ECGM and EMLT models: for these two models only 
one mixing-length parameter is involved and Fc can di- 
rectly be retrieved from a and V. 

KCGM and KMLT models: in the outer region (r < ri) 
as well as in the interior region (t > T2), values of Fc can 
be directly retrieved from a and V. In the transition region 
(ri < T < r2), however, we have to deal with two functions 
for the convective flux: Fi'"* — /i(Vi,ai), the convective 
flux calculated as in the interior, and Fc^^ = hCVujaa), 
the convective flux calculated for the atmosphere. As a 
result of Eq. the convective flux Fc of the model in 
the transition region can be related to Fi'-* and Fi^' as 
follows: 

Fc ^ A(r)Fi^) + (1 - A(t))F« (3) 

where A(r) is - like (3{t) (Eq. ^ - a function of r which 
ensures the continuity of the convective flux. This function 
must decrease with r and must fulfill A(r) = 1 for r < ti 
and A(r) = for t > T2. We point out that both Fc^ and 
Fc'^^ fulfill Fc + Fj-ad — Lq/^ttt^ at any optical depth r. 
Just as for /3(t), the choice for A(r) is rather arbitrary. 
As for the case of V (Eq. 0), we assume A(r) = (3{t) = 

{t2 - r)l {t2 - Ti) for Ti < T < T2. 

Figure |31 shows Fc as a function of t for the KMLT 
and KCGM models. 



Calculation of the driving by the entropy source term 
requires, in addition to the convective flux (Fc), a model 
for the mean anisotropy ($). In the CGM96 model of 
convection the expressions for v'^ and Fc do not depend 
explicitly on the mean anisotropy factor <I>. However, 
CGM96 adopt a model of anisotropy which fixes the ratio 
X = kfi/ky, where and ky are the horizontal and verti- 
cal wavenumbers associated with the eddy of wavenumber 
k (note that fc^ = fc^ -f k^). As a result of that model, for 
the largest scales x = 1/2 while it increases quadratically 
with the total wavenumber k for the smaller scales. Let 
us define $(fc) = v?{k)/u\{k), a k-dependent anisotropy 
factor. For an incompressible fluid - a property assumed 
by the models investigated here (see also CM91) - one can 
show that $(fc) — 1 + x^^. Hence a; = 1/2 implies $ = 3, 
i.e. an isotropic velocity field. As a result of its functional 
dependence on fc, from the larger scales (fc ^ fco) to the 
smaller scales (fc ^ fco) <I>(fc) decreases in the CGM96 
model from ~ 3 towards 1. However, we recall that 
the model of stochastic excitation (MSE) we consider is 
basically isotropic. The anisotropy is taken into account 
only at large scales through Gough's mean anisotropy fac- 
tor ($). Therefore, although CGM96's treatment adopts 
a model where the anisotropy varies with fc, we are left 
with the inconsistency that the turbulent spectrum, F(fc), 
assumed for the MSE (see Sect. I4.1|l is isotropic along 
the turbulent cascade. A possible anisotropy is only be 
taken into account at large scales through $. But as the 
modes are predominantly excited by turbulent eddies with 
fc ~ fco, which carry most of the kinetic energy, this ap- 
proximation is not expected to affect significantly our pre- 
diction. Hence, we assume $ = 3 for the CGM96 models 
in the discussion given below. 

3.2. Convective velocity 

Driving of the oscillation modes by the Reynolds term is 
proportional to pov"^. 

For the CGM96's treatment, is computed according 
to Eq. (88-90) of CGM96. 

The MLT approach provides a theoretical expression 
for w'^ , the mean square of the vertical component of the 
convective velocity. It can be related to the velocity v by 
means of $ (E q.Et. BV's formulation of the MLT assumes 
<b — 2 ("Goueh 19Z7|) . Her e, we compute according to 
Eq. n4Jini of CoxT lilPBSL Chap. 14) and deduce v from 
Eq. 101 with $ = 2. 

Just as the convective flux, the convective velocity v 
can be viewed as a function of V and a, v = /(V, a). 

ECGM and EMLT models: as for the convective flux, 
V can be retrieved from a and V. 

KCGM and KMLT models: in the outer region (i.e. 
T < Ti) as well as in the inner region (r > T2), only one 
a and one V are defined. In those regions v is computed 
as V = /(V, a). However, in the transition region we must 
deal with two different values of a (namely ai and Vi from 
the inner region and aa and Va from the atmosphere) and 
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V is not a linear function of a. Wc thus face the difhcuhy 
to define properly a convective velocity consistent with 
(Eq. 13) in this region. 

We proceed as follows: V(t) and Fc{t) are defined by 
Eq. and Eq. 0, respectively. Then, at fixed V and 
r, wc define an equivalent mixing-length parameter, a*, 
such that Fc — h{V, a*). This parameter is variable with 
depth. We next compute v = /(V(r),Q;*) which hence is 
consistent with Fc of Eq. (O . 

Figure |21 (middle and bottom) shows a* and w as a 
function of depth. For the MLT models, our calculation of 

V assumes <t — 2 which is cons istent with BV's formulation 
of the MLT rsee lCnughllTflTT^ and for the CGM models it 
assumes $ = 3. 
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3.3. Comments 

As shown in Fig. El the EMLT and KMLT models have 
very different convective structures: up until close to the 
top of the superadiabatic region (located at Rq — r ^ 
70 km for those models), the EMLT model results in a 
convective velocity and convective flux smaller than those 
of the KMLT model. This is a consequence of the fact that 



„EMLT 



< a 



KMLT 



(see Table [Hand Table El • On the other 
hand, the EMLT model results in larger v and F^ above 
the top of the superadiabatic region because a 



EMLT 



> 



„KMLT 



a-— = 0.5. 

In contrast with the MLT models, the KCGM and 
ECGM models have rather similar convective structures. 
Indeed, v^cgm ^ecgm (^kcgm F^cgm ^.^gp ) 

have very similar shapes. This in turn is a consequence of 
the fact that the CGM models require values of af^*-'^^ 
and of a^'^'^^ close to that one required for the atmo- 
sphere {Ua = 0.5). 

For the KMLT model, there is a pronounced "kink" at 
the bottom boundary of the transition region {i.e. at t — 
20 or i?Q - r = 100 km), especially for v. For the KCGM 
model the "kink" is much less important. These features 
are directly connected with the angular point observed for 
V in that layer (see Sect. 12.31 and Fig. 

At the bottom of the transition region (i.e. at r = T2), 
for both the KMLT and KCGM models, a*(r) reaches - 
as expected - the "interior" value of the mixing-length 
parameter {i.e. ai), namely: aJ^CGM _ q ^^^^ q,kmlt _ 
2.51. At the top of the transition region {i.e. at r — ri), 
a*(r) reaches for both the KCGM and KMLT models the 
asymptotic value a* ~ 0.50. 

4. Calculation of the solar p mode excitation rates 

4.1. Procedure 

We compute the rate P at which acoustic energy is in- 
jected into solar radial p-modes according to the theoret- 
ical model of stochastic excitation of Paper I and assume 
here - as in Paper II and Paper III - that injection of 
acoustic energy into the modes is isotropic and consider 
only radial p modes. The rate at which a given mode with 
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Fig. 3. Top: Fc is plotted versus Rq — r for the 
KMLT model (dashed curve), EMLT model (dot-dashed 
line), ECGM model (dot-dot-dot-dashed line) and KCGM 
model (solid line). As in Fig.|21the dotted vertical lines de- 
limit the transition region. Middle: Same as the top panel 
for a* (see Sect. 13. 2|) . Bottom: Same as the top panel for 
the root mean square of the convective velocity, v. 



frequency = 27ri^o is excited is then calcu lated with the 
set of Eqs (1-11) of Paper III and Eq. (3) of lSamadi et alJ 
(I2OO4I). 
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The calculation requires the knowledge of three differ- 
ent types of quantities. First of all, quantities which are 
related to the average properties of the medium: 

• the mean density p, as fEa. lA.5(l . the mean square con- 
vective velocity v'^, and the mean square of the entropy 
fluctuations fEa. IA.6jl . They are obtained from the 
equilibrium models as explained in Sect. 13 

Secondly, quantities which are related to the oscillation 
modes: 

• the eigenfunctions and the eigenfrequency v. They 
are cor nputed with the adiaba tic pulsation code 



FILOU llTra,Ti Minh fc T.eonlll 99,^ for each model. 



Finally, quantities which are related to the properties of 
the turbulent flow: 

• the wavenumber (k) dependency of E, i.e. the turbu- 
lent kinetic energy spectrum; 

• the values and depth dependency of fco , the wavenum- 
ber at which the convective energy has its maximum 
and is "injected" into the inertial range of the turbu- 
lent kinetic energy spectrum E; 

• the ^-dependency of Xk, the frequency component of 
the auto-correlation product of the turbulent velocity 
field. 

According to the results in Paper II and Paper III 
obtained on the basis of a 3D numerical solar gran- 
ulation simulation, the /c-dependency of E{k, ly) is ap- 
proximately reproduced by an analytical spectrum called 
'Extended Kolmogoro v Spectrum' (EKS) and defined in 
iMusielak et al.l l)l994l . The z^-dependency of Xk is foimd 
to be better modelled with a Lorentzian function rather 
than by a Gaussian function, which is usually assumed 
for Xk (see Paper III). Within most part of the excitation 
region the spatially averaged anisotropy factor $ is found 
almost constant with a value of $ ~ 2 in agreement with 
BV's value. 

At the top of the superadiabatic region, it was found 
that ko ~ 3.6 Mm~^ and that fco decreases slowly inwards 
with depth (see Paper II). A good approximation for our 
region of interest within the sun is to assume a constant 

4.2. Comparison with the observations 

Results for P are presented in Fig.^ The theoretical esti- 
mates for P are compared with the 'observed' Pobs, calcu- 
lated from iChaplin et al.l l)l998f) 's seismic data according 
to the relation: 



v^^{iy) ^ 27rrMw^(zy), 



(4) 



where i' is the mode frequency, rg is the radius at which 
oscillations are measured. 



is the mode inertia, M = I/i^{rs) is the mode mass, 
and finally P and Vg are the mode line-width and the 
mode surface velocit y, respectively, and are obtained from 
IChapl in et al.l l)l998l) . We point out that, according to the 
definition of Eq. Q), the derived value of Pobs depends on 
the model one considers through M. Indeed, for a given 
mode, M - a priori - differs from one model to another. 
However, the mode masses of the models we consider here 
are very close to each other such that the changes on Pobs 
due to the use of the mode mass of the different models 
are not significant compared to the error bars attached 
to the measurements. For each mode, we can then com- 
pute a unique value for Pobs and compare the values of 
P obtained for the different models to each other and to 
Pobs- We choose to derive M from the radial eigenfunc- 
tions computed for the KCG M model and we ado pt 
rg = Pg-f 200 km consistently with lChanlin et all l|l998>) 's 
observations. 




Fig. 4. The computed solar p mode excitation rates, 
P{v), are plotted versus the frequency for the KMLT 
(dashed fine), EMLT (dot-dashed fine), ECGM (dot-dot- 
dot-dashed line) and KCGM (solid line) models. The filled 
circles represent the 'observed' solar values of P(i') derived 
- according to Eq. Q ~ from the amplitudes and line 
widths of the £ = p-modes measured by iChaplin et alJ 

Ei 



I'M 

1=1 dm^^ 
Jo 



(5) 



As shown in Fig. ^ differences between Pkcgm and 
Pecgm are foun d smaller than the er ror bars associated 
with Pobs from IChaplin et alJ l)l998j) . This is a conse- 
quence of the fact that the ECGM and KCGM models 
present very similar convective structures (see Sect. 13.31 
and Fig. EJ. On the other hand, Pkmlt is found signifi- 
cantly larger than Pemlt as a consequence of the fact that 
the KMLT models result in larger v and Pc values than 
the EMLT ones (see Sect.ESjand Fig.EJl for r > 10 . . . 20. 

Furthermore, Pkcgm and Pecgm are found closer to 
Pobs than Pkmlt- However, above v > 2.5 mHz, differ- 
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ences between Fobs and Pkcgm (or Pecgm) remain impor- 
tant. The origin of this discrepancy is discussed in Sect.El 
For the KCGM and ECGM models we have so far as- 
sumed $ = 3. According to Eq. (jA.6|l . assuming $ = 2 - 
a value which is consistent with results from the 3D so- 
lar simulation - results in a driving by the entropy source 
smaller by a factor ~ 2/3 compared to the case $ = 3. 
This decrease, however, remains small compared to the 
difference in P between the different models and hence 
does not influence our main results. 



1 ^ 



10-6- 



KMLT 
EMLT 
KCGM 
ECGM 



50 100 150 200 250 300 
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Fig. 5. The integrand dP/dm (Eq. (2) of Paper II) is plot- 
ted versus Rq — r for a mode with order n = 2Q in — \ be- 
ing the fundamental mode) for the KMLT model (dashed 
hnc), EMLT model (dot-dashed hne), ECGM model (dot- 
dot-dot-dashed line) and KCGM model (solid line). 



We present in Fig. the integrand dP/dra (Eq. (2) 
of Paper II) corresponding to the excitation rate P of a 
mode of order n = 20. The plot is done for the four solar 
models. For three of the four models the extent of the re- 
gion where most of the excitation takes place is very thin 
(less than 50 km). Obviously, this is the consequence of 
the very shallow extent of the superadiabatic region (see 
Fig- El- This tiny extent of the excitation region strongly 
contrasts with that one found in lStein fc Nordhmdl l|200lh 
and in Paper III. For instance, the latter authors found 
a size of order ~ 500 km for a mode of the same or- 
der. This discrepancy is attributed to the local nature of 
the convective treatments investigated in this work (see 
Sect. El. Note that in Fig. the EMLT model with its 
large value of a for the entire superadiabatic region pre- 
dicts a broader excitation region than the other models. 
This is due to a smaller superadiabatic peak. However, in 
this model the transport of convective heat occurs with 
a smaller convective velocity and hence, as expected, the 
excitation is smaller than for the other models. We note 
here that despite the EMLT model results in a tempera- 
ture structure closer to the numerical simulations (smaller 
superadiabatic peak) , the predicted excitation amplitudes 



are smaller than for the other models investigated here 
and are the most discrepant in comparison with the data. 
This confirms that excitation rates provide a decisive ad- 
ditional test for convection models. 

5. Conclusions 

We have computed the rates P at which acoustic energy 
is injected into the solar radial p modes for several solar 
models. The solar models are computed with two different 
local treatments of convecti on: the cl assical mixing-length 
theory (MLT hereafter) and lCanuto'et al.. (,1996. CGM96 
hereafter) 's formulation. 

For one set of solar models (EMLT and ECGM mod- 
els), the atmosphere is gray and assumes Eddington's ap- 
proximation. The models only assume one mixing-length 
parameter value and reproduce the solar radius at solar 
age but not the Balmer line profiles. For a second set of 
models (KMLT and KCGM models), the atmosphere is 
built using a T{t) law which has been obtained from a 
Kurucz's model atmosphere computed with the same lo- 
cal treatment of convection. The mixing-length parameter 
in the model atmosphere is chosen so as to provide a good 
agreement between synthetic and observed Balmer line 
profiles, while the mixing-length parameter in the interior 
model is calibrated so that the model reproduces the solar 
radius at solar age. 

Both the KMLT and the KCGM models reproduce the 
Balmer line profile and the solar radius and luminosity but 
- as shown in Figs. (2)-(4) - the CGM models model the 
transition between the region of high convective efficiency 
(the interior) and the region of low efficiency (the atmo- 
sphere) in a more realistic way than the MLT models, as 
they reproduce the observed excitation rate P more closely 
and predict a smoother transition region. Furthermore the 
KMLT model requires a change of the mixing-length by 
a factor of five in a layer of ^ 20 km thickness, which is 
significantly less than a pressure scale height 300 km). 
Given the meaning of alpha, this means that the mixing- 
length varies from about 150 km to 750 km in a layer of 
^ 20 km thickness (!), which does not make much sense 
from a physical point of view. On the other hand, the 
KCGM model does not require such large change in a. 

For the MLT treatment, the oscillation excitation 
rates, P, do significantly depend on the properties of the 
atmospheres investigated here. Indeed, differences in P be- 
tween the EMLT model and the KMLT model are found 
to be very large. On the other hand, for the CGM96 treat- 
ment, differences in P between the ECGM and the KCGM 
models are very small compared to the error bars attached 
to the seismic measurements. This result shows that an 
Eddington gray atmosphere can be assumed for the cal- 
culation of P when the CGM96 formulation is adopted. 
This would be particulary convenient in the case of mas- 
sive computations of P for a large set of stellar models. 

For the EMLT and KMLT models, P is signifi- 
cantly underestimated c ompared to th e solar seismic con- 
straints obtained from IChapl m et all lll99al) 's measure- 
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ments. KCGM and ECGM models yield values for P closer 
to the seismic data than the EMLT and KMLT models. 
Contraints on the Balmer line profile and on the so- 
lar radius can be satisfied by the two formulations (MLT 
and CGM96), provided that the mixing- length parame- 
ters ai and are suitably adjusted. Once the above con- 
straints are satisfied, the solar seismic data provide addi- 
tional valuable constraints and according to the present 
investigation (focussed on local approaches) they cleary 
favor the CGM96 treatment. 



Ou r calculations are based here on iGrevesse fc Noelj 
solar abunda nces. There w a s rece n tly a change 
in these values (see lAsplund etal l)2004 l2005|) 'l with 
quite son ie implications for th e internal structur e of th e 
Sun (se e 'Basu & Antia' fiooT); 'Montalba n et all l(2004l) : 
iBahcall et al (2QQa) ; Anti a fc B asu (2005j). Then the im- 
plication for the calculation of the excitation rates (P) 
must be in principe tested consistently by changing the 
solar mixture both in the interior and the atmosphere 
(work in progress). However, we except small changes in 
P. Indeed, as a first step we have calculated a solar model 
having t he low metallicity Z /X = 0.0171 inferred from 
the new lAsplund et al (|20^)'s revised solar abundances 
in which the interior calculation is based on the detailed 
Asplund et al's mixture and associated opacities while the 
atmosphere considers th e low metallicity but keeps the 
iGrevesse fc Noel^ l|l99.'^ 's solar abundances and associ- 
ated opacities. Changes in P smaller than ^ 10% — much 
smaller than the errors bars (^ 20%) associated with the 
current measurements — have been obtained. 

The remaining discrepancy above v > 2.5 mHz be- 
tween computed and observed P (Fig. ^ is attributed 
to the assumption of locality in the present tr e atment 
of convection. As a matter of fact. I Samadi et all l)2003bl 
Paper III) have succeeded in reproducing the seismic con- 
straints much better using constraints from a 3D solar 
granulation simulation. One reason for this improvement 
is that convection is intrinsically a non-local phenomenon. 
In the terminology of classical turbulence modelling, the 
eddies located at different layers contribute to the con- 
vccti ve flux o f a given layer (cf. also the discussion in 
lHoudeklh99(t . Hence, a non-local description of convec- 
tion is expected to predict a more extended superadia- 
batic region. This is suggested, for instance, by the com- 
parison of our present results with that of Paper II. Non- 
local theories - such as tho se by iGous-hl ^ll977^ and by 
[Canuto & Dubovikovl l(l998l) - also support this explana- 
tion: they typically predict a smaller temperature gra- 
dient in the supcradiabatic region than the local theo- 
ries do (see Houdek 1996; Kupka & Montgomerv 2002) 
and thus a more extended supcradiabatic region as well. 
Another property of solar granulation caused by non- 
locality is the observed asymmetry between the areas cov- 
ered by up- and downflows. This allows for a lower root 
mean square velocity while larger velocities (and thus 
more effective mode driving) can be reached in the down- 
flows with their much higher velocities (note that such an 



asymmetry can be accounted for through non-lo cal mod- 
els as proposed by Canuto fc Dubovikovl l)l998l see also 
iKupka fc Montgome ry 2002j). although in a more simpli- 
fied manner). On the other hand, the local models studied 
here cannot account for the different properties of up- and 
downfiow areas at all. Hence, the supcradiabatic region in 
these models is physically different from that one expected 
from non-local models and found in numerical simulations. 
Solar modes above v > 2 mHz, however, are predomi- 
nantly excited in the superadiabatic region (Paper III) . A 
larger amount of acoustic energy is then injected in those 
modes when convection is treated - as is in the 3D sim- 
ulation or on the base of a non-local theory - in a more 
realistic manner than through local theories. 

The results presented here so far only concern our 
Sun. lSamadi~et al. [ 2003^ found that P scales as {L/My 
where s is the slope of the power law and L and M are the 
mass and luminosity of their computed ID stellar models. 
By building a s et of stellar in odels with the MLT and 
another one with lGoughl lll977l) 's non-local formulation of 
convection, the authors found that the slope s is quite sen- 
sitive to the ID treatment of convection. In this respect, 
it would be interesting to compare the influence of t he 
CGM96 formulation or of ICanuto fc Dubovikovl l(l99j^ 's 
non-local convection treatment on the value of s and to 
test whether or not measurements of P can - for stars 
other than the Sun - discriminate the best ID treatment 
of convection. 
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Appendix A: Relation between the entropy 
fluctuations and the convective flux 

As in Paper I, we relate the entropy fluctuations, St, to 
temperature fluctuations, 9, as follows: 



Cp 

-to 



(A.l) 



where Tq is the mean temperature and Cp = {ds/dliiT)p. 
Hence, the mean square of the entropy fluctuations, < 



S( >, can be expressed as 



<0^ > 



(A.2) 



The convective flux is related to 6 and Uz, the vertical 
component of convective velocity, as 



Fc ~ PoCp <Uz9 >, 



(A.3) 
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where pQ is the mean density. We furthermore assume, 
consistently with the adopted quasi-normal approxima- 
tion in Paper I, that < Uz st >^ can be decomposed as 

< 1*2 St >^ = < >< Sj > = < Sj > . (A. 4) 

Finally, one can show that 

as = = PoToFiVad, (A. 5) 

where s is the entropy and p the pressure, Tq is the mean 
temperature, Fi — {dlnp/dln p) ^ is the adiabatic expo- 
nent and Vad = {dhiT/dhip)s is the adiabatic tempera- 
ture gradient. 

The mean square of the entropy fluctuations can then 
be deduced from the set of Eqs. (jA.llA.5|) and Eq. ©: 

where $ is the spatially averaged anisotropy factor, which 
is defined in Eq.®. Driving by the entropy source term is 
then proportional to (as/poTo)^ Fc ~ Paper I 

- and thus scales like the square of Fc- 
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